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I.  INTRODUCTION 

The  research  pursued  under  this  contract  focused  primarily  upon  plasma 
processes  in  and  near  the  polar  cusp.  The  following  sections  will  cover  the 
seven  basic  areas  studied:  Plasma  Injection  and  Transport  in  the  Mid-Altitude 
Polar  Cusp  (Section  II);  Observations  of  Counter-Streaming  Electrons  at  High 
Altitudes  (Section  III);  Observations  of  Upward  Electron  Beams  and  Their 
Relationship  to  Region-1  Birkeland  Currents  (Section  IV);  Observations  of  the 
Electron  Population  Responsible  for  the  6300A  SARarc  emission  (Section  V); 
Polar  Rain  Observations  (Section  VI);  Polar  Wind  Observations  (Section  VII); 
and  Observations  of  Ion  and  Electron  Acceleration  Events  Produced  by  Parallel 
Electric  Fields  (Section  VIII). 

Section  IX  details  the  publications  supported  by  contract,  and  Section  X 
presents  papers  presented  at  scientific  meetings. 

The  primary  observing  platform  for  the  research  reported  here  was 
Dynamics  Explorer  1  (DE-1).  The  DE-1  High  Altitude  Plasma  Instrument  (HAPI) 
consists  of  five  electrostatic  analyzers  mounted  in  a  fan-shaped  angular  array 
at  angles  of  45°,  78°,  90°,  102°,  and  135°  with  respect  to  the  spacecraft  spin 
axis.  Each  analyzer  makes  differential  measurements  of  electrons  and  positive 
ions  over  an  energy/charge  range  of  5  eV/e  to  32  keV/e.  Energy  stepping 
proceeds  at  commandable  rates  of  up  to  64  sec-1,  providing  three-dimensional 
plasma  distribution  functions  at  the  six-second  spin  rate  of  DE-1  (Burch  et 


II.  PLASMA  INJECTION  AND  TRANSPORT  IN  THE  MID-ALTITUDE  POLAR  CUSP 


Ion  fluxes  in  the  cusp  often  show  a  "dispersion"  signature,  in  which  the 
characteristic  ion  energy  fails  with  increasing  latitude  (Shelley  et  al., 

1976;  Reiff  et  al,  1977).  The  dispersion  signature  occurs  because  poleward 
(E  x  B)  convection  acts  as  a  "velocity  filter"  on  the  incoming  plasma.  Par¬ 
ticles  with  large  parallel  velocities  arrive  at  the  earth  sooner  (therefore, 
more  equatorward)  than  particles  witn  lower  parallel  velocities.  Pitch  angle 
effects  on  the  velocity  dispersion  are  not  important  at  low  altitudes  (Reiff 
et  al,  1977).  This  is  not  the  case  in  the  mid-altitude  cusp  (Frank.,  1971; 
Gurnett  and  Frank,  1978),  where  the  transit  time  from  injection  to  observation 
of  a  particle  of  a  given  energy  is  a  strong  function  of  pitch  angle.  This 
fact  enables  us  to  obtain  independent  measurements  of  the  injection  altitude 
for  each  spin  of  the  DE-1  satellite  and  to  determine  local  flow  velocities  as 
low  as  ~10  km/s. 


On  a  typical  polar-cusp  pass,  as  the  satellite  enters  the  cusp,  ion 
fluxes  first  appear  in  higher  (-1-3  keV)  energy  channels  clustered  about  0° 
(downward)  pitch  angles  (a).  Then  the  energy  of  peak  ion  energy  flux  de¬ 
creases  systematically  and  the  distribution  spreads  to  larger  pitch  angles  in 
a  distinctive  'V'  pattern,  in  which  the  peak  energy  flux  occurs  at  higher 
energies  for  larger  pitch  angles.  As  this  pattern  continues  to  evolve,  the 
'V'  signatures  deepen  to  reach  energies  near  200  eV  at  0°  and  spread  through 
all  downwaru  pitch  angles  and  a  wide  range  of  upward  pitch  angles.  Eventual 
iy,  strong  ion  fluxes  are  observed  at  all  pitch  angles  except  for  the  narrow 
loss  cone  (-7°)  around  a  =  180°.  Finally,  a  dropout  develops  at  a~U°  and 
spreads  through  all  downward  pitch  angles  so  that  only  upward  moving  ions, 
which  have  mirrored  below  the  spacecraft,  are  observed. 


We  have  attempted  to  explain  the  ion  energy/pitch-angle  signatures 
observed  by  DE-1  in  the  mid-altitude  cusp  by  application  of  a  simple  plasma- 
injection  model  in  which  cusp  field  lines  are  convected  through  a  spatially 
restricted  region  of  plasma  injection  and  then  into  the  polar  cap.  In  a 
dipole  magnetic  field  the  guiding  center  approximation  is  used  to  compute  an 
E(a,t)  relationship  for  particles  observed  at  a  nominal  geocentric  distance  of 
4  R  and  a  dipole  invariant  latitude  of  78°.  The  relationship  has  the 
folfowing  form: 
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where 


ds  is  arc  length  along  a  dipole  field  line, 
sq,  are  the  observation  and  injection  points, 

M  is  the  particle  mass, 

B(s)  is  the  magnetic  induction  along  the  field  line, 

B  is  the  magnetic  induction  at  the  observation  point, 
o 

a  is  the  observed  pitch  angle, 
o 

t  is  the  transit  time. 

For  a  >  90°,  the  integration  is  carried  down  to  the  mirror  altitude  and  back 

Q 

up  the  field  line  to  the  observation  altitude. 

Since  M  and  t  appear  only  in  the  factor  M _  in  eq.  (1),  the  shape  of  the 

2T2 

log  E  vs.  a0  curve  for  given  injection  and  observation  points  will  be  inde¬ 
pendent  of  transit  time  and  will  be  applicable  to  ions  of  all  masses  as  well 
as  to  electrons.  Computed  log  E  vs.  a0  curves  are  plotted  in  Fig.  1  for 
injection  radii  of  6  Re,  8  R£,  12  Rg,  and  16  Rg,  and  an  observation  radius  of 
4  re*  For  ease  of  comparison  the  four  curves  are  adjusted  to  coincide  at  = 
90°.  Also  plotted  in  Fig.  1  are  the  DE  observations  for  a  single  spin  covering 
the  time  period  from  14:07:30  to  14:07:36  on  September  29,  1981.  The  data 
points  in  Fig.  1  are  strongly  indicative  of  an  injection  distance  in  the 
neighborhood  of  8  Re  (geocentric).  Injection  distances  beyond  12  Re  would 
clearly  not  result  in  the  strong  (E,ao)  dependence  that  is  observed  at  ~4  Re 
by  DE-1. 

As  the  observation  point  moves  to  lower  altitudes,  the  model  predicts 
that  the  log  E  vs.  a0  dependence  will  rapidly  become  weaker,  causing  the 
disappearance  of  any  significant  'V'  shape  at,  for  example,  the  typical  DE-2 
altitude  of  900  km.  This  prediction  is  borne  out  by  the  DE-2  observations  and 
by  the  fact  that  no  other  low-altitude  spacecraft  or  sounding  rockets  have 
observed  the  'V'  shaped  log  E  vs.  a0  relationship. 

Fig.  2  displays  the  computed  time  evolution  of  the  log  E  vs.  a0  relation 
ship  for  H  with  an  injection  point  at  8  R£  and  an  observation  point  at  4  Re> 
All  the  curves  in  Fig.  2  have  the  same  shape,  with  the  vertical  position  of 
each  curve  determined  only  by  the  relationship  vt  =  constant.  The  energy 
scale  in  Fig.  2  is  arbitrarily  truncated  at  just  over  3  keV  and  at  100  eV  to 
coincide  approximately  with  the  high  and  low  energy  limits  of  the  observed  ion 
fluxes.  Fig.  2  then  allows  one  to  trace  the  log  E  vs.  £0  relationship  as  a 
function  of  elapsed  time  after  injection  for  an  input  H  spectrum  that  falls 
off  rapidly  at  energies  above  3  keV  and  below  100  eV^  For  heavier  ions  the 
elapsed  times  in  Fig.  2  would  be  multiplied  by  /M/MH  ;  for  electrons  the  times 
would  be  smaller  by  a  factor  of  43. 
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Figure  1.  Computer  energy/pi tch-angle  relationship  for  particles  injected  at  geocentric 
radial  distance  Rj  on  the  7B°  dipole  magnetic  field  line  and  observed  at  a 
later  time  at  a  distance  of  4  RE  on  the  same  field  line.  The  shaded  region 
is  the  400  km  loss  cone.  Data  points  denote  the  energy  of  peak  ion  energy 
flux  observed  during  one  full  spin  of  DK-1  (dots  for  0°  to  180°  and  tri¬ 
angles  for  180°  to  0°). 
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The  time  evolution  of  the  log  E  vs.  a0  curves  could  result  from  the 
convection  of  cusp  field  lines  through  a  restricted  region  of  plasma  injection 
or  from  a  single  impulsive  injection  event  over  the  entire  cusp.  However,  in 
the  latter  case,  unless  cusp  ion  injection  is  very  infrequent  (several  minutes 
between  injections)  several  of  the  curves  of  Fig.  2  would  exist  simultane¬ 
ously,  and  the  distinct  1 V' -shaped  E,a0  signature  would  not  appear.  Time- 
dependent  injection  no  doubt  occurs  (Carlson  and  Torbert,  1980)  but,  to  be 
consistent  with  DE-1  data,  must  be  confined  to  a  restricted  injection  region 
located  near  the  low-latitude  edge  of  the  cusp. 

The  cusp  velocity  filter  effect  allows  us  to  observe  a  specific  velocity 
at  each  pitch  angle.  The  velocity  selected  by  the  filter  is  gyrotropic,  i.e., 
independent  of  gyrophase  angle,  if  one  is  in  the  rest  frame  of  the  magnetic 
field  line.  If  the  field  line  is  convecting  past  the  spacecraft,  the  convec¬ 
tion  velocity  adds  vectorially  to  the  velocity  selected  by  the  filter  effect. 
Since  each  detector  observes  90°  pitch-angle  particles  twice  per  spacecraft 
spin,  one  can  construct  a  particle  phase-space  distribution  for  the  plane  per¬ 
pendicular  to  the  magnetic  field  direction,  and  have  ten  independent  look 
directions  in  that  plane.  Each  look  direction  will  show  the  velocity  filter 
effect,  with  a  peak  distribution  function  at  a  specific  velocity  within  the 
uncertainty  due  to  the  logarithmic  spacing  of  the  energy  steps.  We  can  then 
fit  an  offset  circle  to  the  velocity  of  the  peak  distribution  function.  The 
magnitude  of  the  offset  gives  us  the  locally  measured  flow  velocity,  and  the 
direction  of  the  offset  yields  the  direction  from  which  the  plasma  is  flowing. 
The  flow  velocities  determined  in  this  way  at  DE-1  are  plotted  along  with  con¬ 
current  DE-2  flow  measurements  (Heelis  et  al,  1981)  in  Fig.  3  for  a  typical 
cusp  pass.  The  agreement  between  the  two  flow  patterns  is  fairly  good,  con¬ 
sidering  the  difficulty  and  uncertainty  in  the  high-altitude  measurement. The 
magnitude  of  flows  in  a  dipole  field  geometry  should  map  upward  along  field 
lines  approximately  as  r3/2.  The  flow  velocities  at  DE-1  should,  therefore, 
be  roughly  six  times  as  high  as  at  DE-2.  Fig.  3  shows  that  the  measured  flow 
velocities  are  more  nearly  15  to  20  times  the  DE-2  flows.  Possible  reasons 
for  the  apparent  enhancement  of  flow  velocities  at  higher  altitudes  will  be 
the  subject  of  further  study. 

The  cusp  plasma  injection  model  presented  above  predicts  that  the  shape 
of  the  log  E  vs.  a0  curves  should  be  the  same  for  ions  and  electrons  if  the 
altitudes  of  injection  and  observation  are  equal.  However,  the  progression  of 
the  curve  down  the+energy  axis  of  Fig.  2  will  proceed  43  times  as  fast  for 
electrons  as  for  H  .  Since  cusp  electron  energies  are  generally  between  10  eV 
and  200  eV,  Fig.  2  shows  that  the  entire  progression  will  take  place  during  a 
time  interval  between  about  5  seconds  and  20  seconds  following  injection.  In 
the  continuous  injection  model  this  prediction  means  that  for  the  'V'-shaped 
signatures  to  be  observed  in  the  electrons,  the  spacecraft  must  cross  a  field 
line  within  about  10  seconds  of  the  time  it  convected  out  of  the  injection 
region.  However,  there  are  times  when  electron  'V's  are  observed.  As 
expected,  the  'V'  signatures  are  observed  in  only  one  or  two  spacecraft  spin 
periods  just  adjacent  to  the  cusp  equatorward  electron  boundary.  The  ion  'V's 
are  observed  to  be  displaced  poleward  from  the  electron  'V's,  again  as 
expected  from  the  simple  model  of  a  spatially  restricted  injection  region  from 
which  cusp  plasma  is  convected  into  the  polar  cap. 


III. 


OBSERVATIONS  OF  COUNTERSTREAMING  ELECTRONS  AT  HIGH  ALTITUDES 


Counterstreaming  electrons  have  been  observed  by  the  High-Altitude  Plasma 
Instrument  (HAPI)  on  the  DE-1  satellite.  Previously,  counter streaming  elec¬ 
trons  have  been  detected  by  the  S3-3  satellite  as  narrowly  collimated  beams  of 
electrons  with  energies  around  one  keV  streaming  both  parallel  and  antipar¬ 
allel  to  the  magnetic  field  direction  (Sharp  et  al.,  1980).  The  S3-3  observa¬ 
tions  were  made  at  altitudes  between  4000  and  8000  km,  while  the  present  work 

pertains  to  observations  of  phenomena  at  higher  altitudes  (~2  R  ). 

E 

Some  features  of  counterstreaming  electrons  suggest  that  the 
mechanism  responsible  for  electron  acceleration  is  not  associated  with  quasi¬ 
static  electric  fields  (Sharp  et  al.,  1980).  For  example,  counterstreaming 
electrons  have  broad  and  unstructured  energy  spectra  that  can  be  fitted  by 
Maxwellian  functions  with  temperatures  varying  from+0.5  to  I  keV. 

Simultaneous  measurements  of  conic  ions  of  H  and  0  in  the  pitch  angle  range 
of  90°  to  125°  dur  ing  counterstreaming  electron  events  also  support  the 
suggestion  that  a  quasi-static  electric  field  is  not  necessarily  involved. 
However,  the  S3-3  satellite  has  occasionally  detected  peaks  in  the  energy 
spectra  of  counterstreaming  electrons. 

Our  survey  of  DE-1  HAPI  data  indicates  that  there  are  two  distinct  types 
of  counterstreaming  electron  events.  The  first  type  of  event  is  fairly  stable 
and  has  broad  and  unstructured  energy  spectra  not  suggestive  of  acceleration 
by  parallel  electric  fields.  The  second  type  is  transient,  and  has  beam-like 
distribution  functions,  which  are  consistent  with  acceleration  by  parallel 
electric  fields.  The  first  type  is  found  to  be  similar  to  those  frequently 
observed  by  the  S3-3  satellite,  while  the  presence  of  the  second  type  is  only 
hinted  at  in  the  S3-3  observations.  These  two  types  of  events  become  quite 
distinguishable  in  the  HAPI  data  because  DE-1  has  a  faster  spinning  rate  and 
the  HAPI  instrument  is  able  to  obtain  the  plasma  distribution  function  at  a 
higher  time  resolution  as  compared  to  the  previous  observations.  We  refer  to 
the  two  types  of  events  as  type  1  and  type  2  counterstreaming  electron  events 
respectively. 

A.  Type  1  Counter streaming  Electrons 

A  typical  example  of  the  electron  distribution  functions  parallel  and 
anti  parallel  to  the  field  direction  measured  during  a  Type  1  counterstream¬ 
ing  event  is  presented  in  Figure  4  (with  X  symbols).  An  electron  distribution 
function  during  one  counterstreaming  electron  event  observed  by  the  S3-3  is 
deduced  from  Sharp  et  al  (1980)  and  presented  here  for  comparison  (with 
triangle  symbols).  Except  for  a  difference  in  the  magnitude,  these  two 
distribution  functions  are  very  similar,  as  both  can  be  approximated  by  a 
superposition  of  two  Maxwellian  functions.  The  temperature  of  the  counter¬ 
streaming  component  is  estimated  to  be  about  200  eV,  while  the  high  energy 
isotropic  component  is  estimated  to  have  a  temperature  of  about  1.5  keV.  It 
appears  that  this  type  of  event  is  related  to  the  counterstreaming  electron 
events  reported  by  Sharp  et  al  (1980). 

The  pitch  angle  range  of  counter  streaming  electrons  observed  is  about 
3bJ,  much  wider  than  the  local  loss  cone  (~id°).  The  counter  streaming 
electron  component  is  found  to  extend  to  larger  pitch  angles  with  smaller 
density;  this  might  be  considered  as  evidence  of  pitch  angle  scattering. 


Figure  4  indicates  that  the  distribution  function  of  counterstreaming 
electrons  is  higher  in  the  downward  direction  <18°<a<30°)  than  in  the  upward 
direction  ( 150° <a<170° ) ,  more  specifically,  about  four  times  higher  at  ener¬ 
gies  less  than  2  keV.  For  this  particular  event,  we  found  the  downward  en¬ 
hancement  to  be  about  2  for  most  times  with  an  occasional  enhancement  in  the 
opposite  direction. 

One  type  of  counterstreaming  electron  frequently  detected  by  the  DE-1 
HAPI  instrument  is  an  enhancement  of  field-aligned  electron  fluxes  in  both 
upward  and  downward  directions.  These  counter streaming  electrons  can  be 
characterized  by  a  Maxwellian  distribution  function  with  a  temperature  of 
about  200  eV.  Their  occurrence  seems  to  be  associated  with  isotropic 
precipitation  of  electrons  with  energies  greater  than  1  keV. 

In  Figure  4  we  have  also  plotted  two  electron  distribution  functions 
of  field-aligned  electron  beams  measured  by  ATS-6  on  the  equatorial  plane 
during  substorm  injection  events  (Lin  et  al.,  1979;  Moore  and  Arnoldy,  1982). 
It  is  interesting  to  notice  that  they  are  characterized  by  two  Maxwellian 
functions  as  well.  The  magnetograms  at  local  midnight  are  not  available  yet 
to  learn  if  type  1  counterstreaming  electrons  occur  during  substorm 
injections,  but  the  Kp  magnetic  activity  index  was  certainly  high  (6-)  for 
this  event.  It  is  possible  that  counterstreaming  electrons  and  substorm 
injected  electron  beams  are  related. 

A  simple  interpretation  of  type  1  counterstreaming  electrons  is  that 
these  electrons  are  produced  at  the  topside  ionosphere  as  thermal  electrons 
heated  to  higher  temperature  through  wave-particle  interactions.  The  heated 
electrons  appear  as  streaming  particles  as  they  move  up  along  field  lines. 

In  order  to  explain  the  counterstreaming  feature,  the  heating  should  take 
place  continuously  for  long  duration  above  both  ionospheres.  However,  this 
interpretation  does  not  explain  the  events  described  by  Sharp  et  al  (1980) 
which  were  much  narrower  than  the  loss  cone. 

Another  possible  interpretation  is  that  type  1  counterstreaming 
electrons  are  secondary  electrons  produced  by  collisions  of  precipitating 
electrons  with  the  atmosphere,  as  suggested  for  substorm-injected  electron 
beams  by  Moore  and  Arnoldy  (1982).  Wave  particle  interactions  then  cause  the 
electrons  to  be  scattered  outside  the  loss  cone.  The  diffusion  coefficient  D 
can  be  estimated  according  to  (Aa)2/2At,  where  Aa  is  the  average  pitch  angle 
change  and  At  is  the  interaction  time.  We  estimate  Aa  to  be  about  12°  and  At 
to  be  the  elapsed  time  for  electrons  of  500  eV  to  travel  from  the  topside 

ionospheres  to  the  satellite.  We  find  D  to  be  about  9  x  10-3  sec-1  for 

electrons  originating  from  the  northern  hemisphere  and  3  x  10-3  sec-1  for 

electrons  originating  from  the  southern  hemisphere.  The  estimated  D,  which 

represents  an  upper  limit,  agrees  in  order  of  magnitude  with  the  diffusion 
rate  estimated  for  substorm-injected  electron  beams  (Lin  et  al.,  1979;  Moore 
and  Arnoldy,  1982). 

We  notice  that  the  electron  distribution  function  observed  by  the  S3- 
3  is  about  one  order  of  magnitude  larger  than  the  DE-1  measurement.  This 
difference  could  be  attributed  to  different  events  or  instrumental  effects. 
This  difference  might  also  be  interpreted  as  an  altitude  dependence  since  the 
S3-3  altitude  was  much  lower  than  the  DE-1  altitude. 
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B.  Type  2  Counters tr earning  Electrons 

Contours  of  the  logarithm  of  the  electron  distribution  function  are 
plotted  in  Figure  5  for  three  6-second  time  intervals  during  a  type  2  event: 
(a)  before,  (b)  during  and  (c)  after  the  counterstreaming  electron  event.  In 
the  figure,  solid  lines  of  circles  and  hyperbolas  are  theoretical  predictions 
of  the  demarcation  boundary  between  ionospheric  and  magnetospheric  populations 
when  a  quasi-static  electric  field  is  assumed  above  and  below  the  satellite 
(Chiu  and  Schulz,  1978;  Mizera  et  al.,  1981).  Following  Mizera  et  al.  (1981), 
we  infer  the  potential  drop  below  the  satellite  to  be  the  ion  beam  energy  and 
the  potential  drop  above  the  satellite  to  be  the  electron  beam  energy.  The 
electron  and  ion  beam  energies  are  determined  from  the  plot  of  the 
distribution  function  along  the  axis.  The  computed  demarcation  curves 
appear  to  characterize  reasonably  the  phase  space  density  contours  before  and 
after,  but  not  during  the  event.  More  specifically,  the  hyperbolas  in  Figures 
5a  and  5c  mark  well  the  electron  loss  cone  and  the  low  energy  core  of 
presumably  ionospheric  electrons.  Figures  5a  and  5c  indicate  a  slight 
enhancement  of  phase  space  density  in  the  forbidden  region  between  the  circle 
and  the  hyperbola.  These  features  resemble  the  S3-3  observations  of  electron 
phase  space  density  contours  during  electrostatic  shocks  (e.g.  see  a  review 
by  Fennell  et  al.,  1981).  In  contrast,  the  hyperbola  in  Figure  5b  lies  well 
inside  the  region  of  the  low  energy  core.  Furthermore,  Figure  5b  displays  a 
distinct  counterstreaming  feature,  an  electron  beam  at  an  energy  of  416  eV 
( v ||  -  12000  km/sec)  along  the  +V||  axis  and  a  weak  electron  beam  at  an  energy 
of  about  100  eV  (v  =  -5920  km/sec)  along  the  -v.j  axis.  The  phase  space 
density  contours  snow  no  loss  cone  feature  and  are  very  elongated  along  the  v 
axis  up  to  about  v^  =  15000  km/sec  (or  800  eV).  Figure  5b  therefore  cannot  be 
explained  by  a  simple  potential  drop  that  increases  monotonically  from  the 
magnetosphere  to  the  ionosphere. 

The  type  2  counterstreaming  electron  event  is  distinguished  from  type 
1  events  in  several  respects.  First  of  all,  the  electron  distribution 
function  has  beam-like  structure.  Second,  the  electron  energy  flux  lacks  a 
hot  component  (>1  keV) .  Finally,  the  ion  energy  fluxes  observed  during  the 
type  2  event  have  pitch  angles  near  175*,  much  closer  to  the  magnetic  field 
direction  than  those  of  the  conic  ions  detected  in  the  type  1  events.  The 
type  2  event  may  have  been  observed  by  the  S3-3  satellite,  since  S3-3 
occasionally  detects  energy  spectra  with  beams  during  counterstreaming 
electron  events  (Sharp  et  al.,  1980). 


One  unusual  feature  of  type  2  counterstreaming  electrons  is  the 
elongation  of  the  phase  density  contours  along  the  v  axis  (see  Fig.  5b).  The 
elongated  contours  can  be  interpreted  as  electron  beams  reflected  by  a 
potential  barrier.  The  HAPI  measurement  therefore  suggests  an  upward  electric 
field  above  the  satellite  and  a  downward  electric  field  below  the  satellite. 

To  produce  such  an  oppositely-directed  electric  field,  we  suggest  V-shaped 
potential  contours  imbedded  with  a  potential  island  (valley).  A  schematic 
illustration  is  given  in  Fig.  6.  This  potential  profile  could  explain  the 
decrease  of  ion  beam  energy  if  the  satellite  passed  through  the  transition 
region  between  the  po  tential  island  and  the  V-shaped  potential  as  shown  in 
Fig.  6. 
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Figure  5. 


Electron  phase  space  density  contours  (a)  before,  (b)  during  and 
(c)  after  a  type  2  counterstreaming  electron  event. 


Sharp  et  al  (I960)  suggested  another  interpretation  of  counter¬ 
streaming  electrons  in  terms  of  a  double  layer  with  a  potential  drop  of  about 
1  keV  flickering  its  polarity  about  every  1  millisecond.  To  account  foe  the 
HAPI  observation,  we  estimate  the  flickering  period  to  be  at  least  50  milli¬ 
seconds  since  the  counterstreaming  portion  of  the  electron  spectra  was 
detected  successively  for  three  energy  steps.  Low  altitude  satellite 
observations  of  auroral  particle  fluxes  have  not  yet  suggested  flickering 
double  layers  of  such  long  periods. 

The  observed  sequence  of  the  type  2  event  presented  in  Figure  5  could 
be  interpreted  either  as  a  spacecraft  horizontally  crossing  a  V-shaped 
potential  structure  or  as  a  vertical  motion  of  the  structure  passing  the 
spacecraft  along  magnetic  field  lines  toward  the  ionosphere.  If  DE-1  were 
detecting  a  spatial  structure,  the  transverse  width  of  the  potential  island 
could  be  as  narrow  as  16  km  along  the  satellite  path  which  corresponds  to  a 
North-South  distance  of  about  1  km  at  the  ionosphere.  On  the  other  hand,  if 
DE-1  were  detecting  temporal  variations,  the  potential  island  would  be  very 
transient,  lasting  only  about  one  satellite  spin  period. 

Potential  contours  with  an  island  structure  have  recently  been  found 
in  a  computer  simulation  of  two-dimensional  double  layers  in  which  the  effects 
of  magnetic  mirror  force,  trapped  electrons  and  background  plasmas  are 
included  (Wagner,  1980;  Wagner  et  al.,  1981).  Recent  laboratory  experiments 
on  double  layers  also  indicate  a  similar  structure  (Iizuka  et  al.,  1982).  The 
condition  for  producing  such  a  potential  structure  is  not  clear,  although 
computer  simulations  by  Wagner  (1980)  suggest  strong  wave  trapping  due  to  the 
Buneman  instability.  We  notice  that  the  potential  profile  in  Figure  6  is 
very  similar  to  the  magnetic  field  and  electric  potential  overshoots  of  the 
earth's  bow  shock  which  arise  from  reflections  of  a  portion  of  the  solar  wind 
ion  population  (Russell  and  Greenstadt,  1979;  Leroy  et  al.,  1981)  .  This 
analogy  might  be  applied  to  the  present  case  if  we  speculate  that  the 
potential  well  in  Figure  6  is  caused  by  reflection  of  ionospheric  electrons. 

Finally,  we  note  that  the  width  of  the  loss  cone  shown  in  Figures  5a 
and  5c  is  energy  dependent  (about  33°  at  738  eV  and  about  40°  at  416  eV). 

These  loss  cone  widths  are  found  to  agree  with  the  presence  of  a  potential 
drop  below  the  spacecraft  according  to  the  formula  (Sharp  et  al.,  1979) 

1  1/2 

a  =  sin-  (B  (E  +  eV)/B.  E)  (2) 

s  1 

where  E  is  the  energy,  eV  is  the  potential,  B  and  B.  are  the  magnitudes  of 
the  magnetic  field  at  the  locations  of  the  satel liteXand  the  top  of  the  iono¬ 
sphere  respectively. 

C.  Correlative  Observations  of  Counter streaming  Electrons 

1 .  Plasma  Measurement 

The  plasma  data  obtained  by  the  High  Altitude  Plasma  Instrument 
(HAPI)  l Burch  et  al^,  1981]  on  DE-1  are  shown  in  Figures  12-16  for  18  September 
1981  (day  261).  Figure  s  a  spectrogram  of  electron  energy  fluxes  for  a  26- 
minute  time  interval  fro  zl28  to  2154  UT.  The  top  panel  shows  electron 
energy  fluxes  at  0°  pitch  angle  (downward),  and  the  bottom  spectrogram  shows 
electron  energy  fluxes  at  90°  pitch  angle.  The  vertical  scale  is  energy  in  a 
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logarithmic  scale,  and  the  horizontal  scale  is  UT.  The  energy  fluxes 
measured  at  both  pitch  angles  were  enhanced  after  2133  UT,  indicating  that  DE- 
1  encountered  a  plasma  injection  boundary.  At  this  time,  the  spacecraft  was 
at  67s  invariant  latitude,  2400  (MLT)  and  14,000  km  altitude,  and  moving 
eguatorward  near  the  midnight  meridian  plane.  A  second  enhancement  of  energy 
flux  occurred  after  2146  UT.  The  energy  flux  remained  at  high  levels  until 
the  data  gap  at  2154  UT. 

This  summary  spectrogram  shows  two  distinct  populations  of 
electron  fluxes.  At  energies  around  several  keV  are  the  energetic 
precipitating  electron  fluxes  that  are  believed  to  originate  in  the  plasma 
sheet.  At  lower  energies,  from  10  to  several  hundred  eV,  enhancements  of 
energy  fluxes  were  detected  sporadically  around  2135,  2139,  and  2148  UT. 

These  low  energy  electron  fluxes  are  enhanced  in  the  field-aligned  direction 
because  the  observed  energy  flux  was  larger  at  0s  pitch  angle  (top  panel)  than 
at  90s  pitch  angle  (bottom  panel).  These  enhanced  electron  fluxes  in  the 
field-aligned  direction  are  similar  to  those  observed  at  synchronous  orbit 
during  plasma  injection  events  [Lin  et  al . ,  1979;  Moore  and  Arnoldv.  1982]. 

In  Figure  8,  we  show  the  electron  density  and  energy  flux 
integrated  from  18  eV  to  10  keV  over  the  whole  pitch  angle  range.  Before  the 
plasma  injection  around  2133  UT,  the  number  density  was  about  1  cm-3  and  the 
integrated  energy  flux  was  about  1  ergs/cm2-s.  After  the  injection,  the 
number  density  increased  to  about  10  cm-3  and  the  integrated  energy  flux 
increased  to  10  ergs/cm2 -s.  Electrons  with  energy  less  than  1  kev  contributed 
about  10%  of  the  integrated  energy  flux  and  about  70%  of  the  number  density. 
Around  ^148  UT,  the  integrated  energy  flux  reached  a  peak  value  exceeding  100 
ergs/cm  -s.  We  found  this  plasma  injection  event  to  be  the  most  intense  in 
the  HAPI  data  that  have  been  processed  so  far.  The  three-hour  Kp  index  during 
this  time  was  6-. 

The  electron  and  ion  energy  fluxes  for  the  full  pitch  angle  range 
( 0* -  180*)  are  presented  in  Figure  9  for  one  counterstreaming  electron  event 
during  2147-2150  UT.  The  DE-1  altitude  in  this  time  interval  was  about  11400 
km.  The  energy  fluxes  shown  are  spin-modulated  with  a  period  of  6  seconds. 

The  electron  spectrogram  (top  panel)  shows  bright  vertical  stripes  below  1 
keV,  indicating  an  enhancement  of  energy  flux  in  that  energy  range.  The 
center  panel,  which  displays  the  detector  pitch  angle,  shows  that  the  detector 
measured  enhanced  energy  fluxes  only  near  0*  and  180*.  Hence  the  measured  low 
energy  (<1  keV)  electrons  were  counterstreaming  along  field  lines.  This 
counter  streaming  electron  event  is  classified  as  type  1,  because  energetic 
electron  fluxes  above  1  keV  were  simultaneously  detected. 

Ion  fluxes  from  approximately  10  to  100  eV  were  detected  at  the 
same  time  that  counters treaming  electrons  were  observed  (see  the  bottom  panel 
of  Figure  9).  In  conjunction  with  the  enhancement  of  energetic  electron 
fluxes,  the  energy  of  conic  ions  increased  from  about  10  eV  to  about  100  eV. 
These  ion  fluxes  are  identified  as  conic  ions  because  their  pitch  angles  lay 
between  110*  and  160*.  The  absence  of  field-aligned  ion  fluxes  suggests  that 
the  observed  conic  ions  are  not  accelerated  by  parallel  electric  fields  below 
the  spacecraft. 


Energy-time  spectrograms  of  electron  energy  fluxes  at  0°  pitch  angle  (top)  and  90°  pitch  angle  (bottom) 
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In  Figure  10,  we  present  contours  of  constant  distribution  • 

function  for  a  12-second  interval  (2147:40-2147:52  UT)  during  which  the 

counterstreaming  electron  flux  was  intense.  The  contour  plots  are  prepared  \ 

from  a  two  full  spin  period  of  data.  The  locus  of  the  sampled  data  in  the  | 

phase  space  starts  at  the  upper  right  corner  and  is  shown  as  dot.  The  | 

contours  at  energies  less  than  1  keV  are  elongated  along  the  axis.  An  J 

island  contour  at  V..  K  9250  km/s  indicates  a  downward-moving  beam  at  an  energy  , 

of  about  250  eV  (positive  velocity  downward).  As  in  this  case,  electron  beams  ( 

are  sometimes  found  when  counterstreaming  electron  fluxes  are  intense.  But  ! 

when  counterstreaming  electron  fluxes  are  weak,  electron  beams  are  generally  j 

not  observed.  Because  the  plasma  data  are  dominated  by  spacecraft 
photoelectrons  below  15  eV,  we  cannot  determine  whether  counterstreaming 
electrons  have  a  beam  energy  less  than  15  eV. 

In  Figure  11  we  show  two  consecutive  contour  plots  of  constant 
distribution  function  that  might  be  considered  as  representative  for  j 

counter streaming  electron  events.  The  two  contour  plots  are  in  time  sequence, 
and  each  plot  is  composed  of  the  plasma  data  during  a  two-spin  period.  The 
consistent  feature  is  that  the  contours  are  elongated  along  the  '•  axis. 

Figure  11  differs  from  Figure  10  in  that  the  contours  do  not  show  electron 
beam.  In  general,  electron  beams  are  found  only  when  counterstreaming 
electron  fluxes  are  intense.  Figure  11  also  shows  that  the  contours  are  less 
elongated  than  the  contours  shown  in  Figure  10  which  covers  a  time  period 
about  10  seconds  earlier.  This  result  suggests  that  ♦'he  angular  width  of 
field-aligned  electron  fluxes  can  vary  in  time.  From  the  intersections  of 
contours  of  constant  distribution  with  the  two  velocity  axes,  we  note  that  the 
contours  in  Figure  16  are  slightly  more  elongated  in  the  direction  of  positive 
velocity  ( toward  the  ionosphere).  From  a  survey  of  contour  plots,  we  conclude 
that  the  counterstreaming  electron  fluxes  can  be  slightly  enhanced  in  either 
upward  or  downward  direction. 

The  magnetic  field  measurements  during  the  counterstreaming 
electron  events  on  day  261,  1981,  are  shown  in  Figure  12.  The  magnetic  field 
instrument  was  described  in  detail  by  Farthing  et  al . ,  [1981].  The  earth's 
internal  field  has  been  subtracted  from  the  observed  magnetic  field.  The 
residual  field  components  were  then  transformed  into  DBR ,  DBTHETA,  and  DBPHI 
in  geomagnetic  dipole  coordinates,  where  DBR  is  on  the  radial  direction. 

DBTHETA  is  perpendicular  to  the  radius  of  the  meridian  plane,  and  DBPHI 
points  toward  the  west.  The  magnetic  field  data  plotted  is  1/2-second 
averaged.  The  magnetic  field  disturbance,  which  was  mainly  in  the  DBPHI 
component,  began  around  2132  UT. 

2 .  Magnetic  Field  Measurements 

Figure  7  indicates  tnree  counterstreuming  electron  events 
occurring  in  the  intervals  21j4-213o,  21JB-2140  ana  UT.  Magnetic 

field  perturbations  in  the  azimuthal  direction  were  associated  with  these 
three  events.  During  the  counterstreaming  electron  event  (2147-zl5u  UT), 

DBPHI  decreased  by  200  nT  in  20  seconds  starting  from  2147:15  UT  (Figure  12). 

This  steep  decrease  of  DBPHI  represents  field-aligned  currents  flowing  into 
the  ionosphere.  From  2148  to  2149  UT,  DBPHI  slowly  recovered  to  its 
unperturbed  value.  The  energy  flux  of  counterstreaming  electrons  was  most 
intense  (2147:30  -2148  UT )  wnen  DBPHI  was  at  its  minimum,  rather  than  at  its 
steep  gradient. 


(KM/S) 


DE-1  HflPI  ELECTRONS 
81261  (18  SEP.  81) 
21:47:40.0  TO  21:47:52.2  U.T. 


DISTRIBUTION  CONTOURS  LOG  F(SEC*/KH*I 


ILRT  101 

6  2.2 

H  IN  EV  HR*  EV 

RV.  Ev. 

AHPS/R* 

ML  T  mft> 

2 3.  7 

AIT  IKHI 

11620 

6.  9900. 

H . I2E*02 

3.66E-07 

GIRT  (Ol 

HO.  S 

GlONG  (O' 

16.3 

FLU*  l /CN*$) 

DEN  1 ✓ CM*  j 

ERGS/CH»S 

0  -  90 

0  -  90 

0  -  90 

ILRT  (Ol 

62.  1 

l.$3E*10 

l . 98E*0l 

3. 67E  *01 

HI  T (HR) 

23.  7 

RLT  IKHI 

1  157* 

GIRT  (01 

HO .  3 

VRRR  IKN/S* 

PC  1HHOS1 

HC (NHOS) 

GLONG IDi 

16.2 

5.69E*01 

N/R 

N/R 

FIGURE  10. 


Phase-space  density  contours  of  the  logarithm  of  the  electron 
distribution  function.  A  downward  electron  beam  is  displayed 
as  a  contour  island  centered  on  the  positive  Vm  axis. 
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FIGURE  11 


Phase-space  contours  of  constant  distribution  function  during  a 
type-1  counterstreaning  electron  event. 
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PIGURE  12.  Three  components  of  magnetic  fields  for  the  DE-1  pass  shown  in  Figure  1 
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3 .  STARE  Measurements 


On  this  pass,  DE-1  was  magnetically  conjugate  with  STARE. 

STARE's  geographic  longitude  is  16°  east,  very  close  to  the  DE-1  orbital 
plane.  The  STARE  radars  recorded  the  intensity  and  Doppler  shifts  of  radio 
waves  scattered  by  E-layer  electron  density  irregularities  [ Nielsen  and 
Greenwald,  1979],  By  combining  the  Doppler  shifts  observed  from  two  separate 
radars,  one  then  estimates  the  electron  E  X  B  drift  velocity  of  the  E-layer 
irregularities . 

Figure  13  illustrates  the  latitudinal  profile  of  electron  drift 
velocities  observed  by  the  STARE  system  from  2120  to  2143  UT.  The  electron 
drift  velocities  are  plotted  in  a  format  that  has  geomagnetic  latitude  as  the 
vertical  coordinate  and  UT  as  the  horizontal  coordinate.  Each  vector 
represents  a  one-minute  average  of  radar  measurements.  The  electron  drift 
velocity  started  to  increase  at  about  2130  UT  in  the  region  equatorward  of  68° 
geomagnetic  latitude,  about  three  minutes  before  the  plasma  injection  boundary 
was  detected  by  DE-1.  Because  DE-1  was  slightly  southward  of  STARE,  this 
observation  can  be  interpreted  by  an  equatorward  motion  of  the  plasma 
injection  boundary.  Since  most  drift  vectors  were  directed  eastward,  the 
electric  field  (E  =  -V  X  B)  was  directed  equatorward.  This  suggests  that  a 
Pedersen  current  was  flowing  in  the  equatorward  direction.  The  magnitude  of 
the  Pedersen  current  increased  with  latitude  because  the  drift  velocity 
increased  with  latitude.  The  increase  of  Pedersen  current  with  latitude 
suggests  an  upward  field-aligned  current  flowing  in  the  observation  area  of 
STARE  after  2130  UT. 

The  field-aligned  current  density  is  estimated  from  STARE  data. 
Assuming  a  homogeneous  Pedersen  conductivity  of  1  mho,  we  estimate  that  the 
upward  field-  aligned  current  density  at  2133  UT  was  about  0.1  uA/m2.  At 
2148  UT,  the  field-aligned  current  was  upward  but  flanked  by  downward 
currents.  The  magnitude  of  field-aligned  current  density  is  estimated  to  be 
about  0.4  yA/rn2  for  a  homogeneous  Pedersen  conductivity  of  1  mho.  Due  to 
uncertainties  on  the  electric  field  magnitude  derived  from  STARE  data,  the 
current  density  estimates  are  very  rough.  However,  the  directions  of  field- 
aligned  current  deduced  from  STARE  are  usually  reliable. 

4.  Ground  Magnetograms 

The  H-components  recorded  by  five  ground  magnetograms  in  the 
northern  auroral  zone  are  shown  in  Figure  14.  These  magnetograms  indicate 
that  a  magnetic  substorm  started  around  1940  UT  at  the  stations  in 
Scandinavia.  Near  21  hr  UT,  all  magnetograms  recorded  active  magnetic 
disturbances.  According  to  the  Abisko  and  Tromso  stations  at  about  2°  east  of 
the  DE-1  orbital  plane,  the  H  component  began  to  recover  after  2130  UT  and 
returned  to  its  pre-substorm  value  after  22  hr  UT  at  Tromso  and  after  23  hr  UT 
at  Abisko,  We  conclude  that  the  plasma  injection  event  observed  by  DE-1 
occurred  during  the  recovery  phase  of  a  magnetic  substorm. 

5 .  Computation  of  Current  Density 

We  compute  field-aligned  current  density  from  the  magnetic  field 
data  and  from  the  plasma  data.  Assuming  a  series  of  field-aligned  infinite 
current  layers  that  are  locally  normal  to  the  dipolar  meridian  plane,  the 
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FIGURE  13 


STARE  latitudinal  profile  of  electron  drift  velocities  versus 
tine  on  day  261,  1981. 


azimuthal  (east-west)  component  of  the  magnetic  field  DBPHI  is  differentiated 
to  yield  the  current  density.  The  differentiation  of  tne  magnetic  field  is 
made  with  respect  to  the  distance  normal  to  the  current  sheets.  From  the 
measured  electron  distribution  f(v,a,<p),  the  field-aligned  current  density  is 
calculated  by 
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where  a  is  pitch  angle, 


e  is  electron  charge. 


and  V  are  the  velocities 
1 


parallel  and  perpendicular  to  magnetic  field  B,  V  .  and  V  are  the 

min  max 

velocity  limits  on  the  measurement,  and  d  is  the  azimuthal  angle  about  B.  We 
set  V  equal  to  the  velocity  of  an  18  eV  electron  to  avoid  the  spacecraft 

photoelectrons,  and  set  V  ^  equal  to  the  velocity  of  a  10  keV  electron. 

The  time  interval  to  acquire  the  electron  distribution  function  over  the 
whole  pitch  angle  range  is  6  seconds  for  the  normal  mode  operation  of  the 
plasma  instrument. 


Figure  15  illustrates  the  current  density  computed  for  the 
interval  corresponding  to  Figure  9.  The  current  density  computed  from  plasma 
data  is  given  in  the  top  panel  and  the  current  density  computed  from  magnetic 
data  is  given  in  the  center  panel.  The  azimuthal  component  of  the  magnetic 
field  is  also  presented  at  the  bottom  panel  for  comparison.  The  plasma  data 
suggest  a  downward  current  density  with  a  peak  value  of  1.5  pA/m2  occurring  at 
2147:15  UT.  The  peak  current  density  derived  from  the  magnetic  field  data  is 
about  2  pA/m2.  Integrated  over  the  current  layer,  the  total  current  densities 
deduced  from  both  methods  are  in  good  agreement  and  are  approximately  0.1  A/m 
within  10%. 


The  magnetic  field  data  indicate  upward  currents  on  both  sides  of 
the  downward  current  with  a  peak  current  density  around  0.5  pA/m2 .  This 
current  density  is  consistent  with  the  estimate  based  on  the  STARE  data,  but 
the  plasma  data  suggest  a  large  upward  current  around  2148  UT  with  a  peak 
density  of  about  5  uA/ra2 ,  about  ten  times  larger  than  current  density  computed 
from  the  magnetic  data.  The  disagreement  could  mean  that  the  upward  current 
was  a  filament  or  that  the  upward  plasma  current  was  partially  neutralized  by 
cold  electrons  moving  in  the  upward  direction. 

The  contribution  of  field-aligned  current  density  by  counter¬ 
streaming  electrons  was  investigated  statistically.  We  surveyed  the  available 
DE-1  plasma  data  that  have  been  processed  so  far,  and  found  a  total  of  18 
counter  streaming  electron  events.  These  events  were  observed  in  the  evening 
auroral  zone,  where  their  magnetic  local  times  were  between  20  and  24  hr,  and 
their  invariant  latitudes  varied  from  60°  to  72*.  These  events  occurred  at 
altitudes  between  9000  km  and  14500  km. 

We  computed  differential  current  density  below  235  eV, 

I  ( <2 J5  eV),  and  the  total  current  density,  I  from  the  plasma  data.  The 

tOta  X 


stepped  energy  of  235  eV  was  chosen  as  the  threshold  because  the  energy  of 
counterstreaming  electrons  is  found  to  be  less  than  250  eV  for  most  events. 
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PIGURE  14. 


Ground  magnetograms  in  the  northern  auroral  zone  on 
18  September  (day  261),  1981. 


The  survey  of  S3-3  data  also  indicates  that  70%  of  electron  beams  were 

detected  below  240  eV  [Collin  et  al.,  1982].  The  largest  current  density 

found  was  5  uA/m2  on  day  81261.  Most  counterstreaming  electron  events  had  a 

total  current  density  less  than  2  pA/ra2.  The  total  current  density  was 

smaller  than  0.25  yA/m2  for  only  two  events.  For  those  events  with  I 

1  total 

greater  than  0.25  uA/m2 ,  we  plot  the  ratio  I(<235  eV)/I  ,  versus  I  ,  in 

total  total 

Figures  16  and  17.  For  the  downward  current  (open  circle),  I (<235  eV)  was 
greater  than  40%  of  the  total  current  density  (Figure  16).  For  the  upward 
current  (solid  circle),  I(<235  eV)  was  smaller  than  50%  of  the  total  current 
density  (Figure  17). 

6.  Energy  and  Pitch  Angle  Distributions 

In  this  section,  we  examine  the  properties  of  electron 
distribution  functions  during  type  1  counterstreaming  electron  events.  We 
present  line  plots  of  distribution  functions  with  small  pitch  angles  (<25° )  in 
the  energy  range  of  20-1000  eV  during  the  interval  from  2147:30  to  2149  UT 
(Figure  18).  To  the  first  order  approximation,  these  distribution  functions 
can  be  represented  by  Maxwellian  functions,  but  weak  electron  beams  with 
energy  less  than  200  eV  were  also  detected,  most  distinctly  in  Figure  23b,  c, 
and  g.  These  electron  beams  moved  in  the  downward  direction  and  were  observed 
in  the  region  of  upward  field-aligned  currents  (Figure  20).  The  distribution 
functions  for  electrons  with  pitch  angles  greater  than  155°  for  the  interval 
2147  -2149  UT  are  shown  in  Figure  19.  Only  two  out  of  the  eight  distribution 
functions  shown  contain  weak  electron  beams  with  energies  around  300  eV 
(Figure  19  and  c).  These  upward  electron  beams  appeared  to  be  insignificant 
in  determining  the  direction  of  field-aligned  current  because  they  were 
detected  in  the  region  of  upward  field-aligned  currents  (Figure  15). 

D.  Conclusions 

Type  1  counterstreaming  electrons  are  apparently  the  type  of  electron 
event  first  reported  by  Sharp  et  al.  (1980),  and  may  be  related  to  the  field- 
aligned  electron  beams  observed  by  the  ATS-6  satellite  at  synchronous  orbit 
during  substorm  injections  (Lin  et  al.,  1979;  Moore  and  Arnoldy,  1982).  Type 
2  counterstreaming  electrons  are  a  new  and  unusual  phenomenon  that  DE-1  has 
detected  within  acceleration  regions.  The  main  feature  distinguishing  these 
two  types  of  event  is  that  the  electron  distribution  function  of  type  1  is 
composed  of  two  Maxwellian  distributions  while  the  distribution  function  of 
type  2  is  characterized  by  electron  beams. 

The  presence  of  two  distinct  types  of  counterstreaming  electron 
events  indicates  the  possibility  of  two  distinct  mechanisms  of  electron 
acceleration  operative  at  high  altitudes  in  the  auroral  zone.  Type  1  events 
appear  to  involve  wave-particle  interactions  whereas  type  2  events  seem  to 
result  from  direct  acceleration  by  oppositely-directed  electric  fields 
pointing  toward  the  satellite  along  magnetic  field  lines. 
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Figure  16.  Ratio  of  current  densities  I ( <235eV)/Itotai  versus  Itotal  for 

downward  currents.  The  current  density  I(<235eV)  is  integrated 
from  the  plasma  data  from  18  to  10,000eV. 
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IV.  UPWARD  ELECTRON  BEAMS 


Introduction  and  Observations 


The  carriers  of  the  large-scale  Birkeland  currents  linking  the 
ionosphere  and  the  magnetosphere  have  typically  been  thought  to  be  upward- 
moving  thermal  electrons  in  the  downward-current  regions  and  precipitating 
auroral  electrons  in  the  upward-current  regions  (Anderson  and  Vondrak,  1975; 
Klumpar ,  1979).  However,  the  problem  is  complex,  and  both  populations  can 
coexist  in  regions  of  upward  or  downward  net  current.  For  example,  in  the 
duskside  plasma  sheet  where  the  current  is  typically  downward,  the  upward  flux 
of  cold  ionospheric  electrons  must  be  larger  than  that  inferred  from  in  situ 
magnetometer  measurements,  since  it  must  also  balance  the  upward  current 
caused  by  the  precipitation  of  energetic  electrons.  Similarly,  intense 
auroral  electron  precipitation  is  often  accompanied  by  upward  fluxes  of  near- 
thermal  and  low-energy  (50  eV)  electrons  whose  current  can  balance  roughly 
half  of  the  primary  upward  current  (Maier  et  al^,  1980).  These  low-energy 
electrons  are  probably  secondaries  produced  by  the  auroral  beam  itself 
(Evans,  1974;  Pulliam  et  al,  1980). 

Narrowly-collimated  electron  beams  also  can  contribute  significantly 
to  the  total  current  balance.  Unlike  the  substantially  isotropic  auroral 
electron  precipitation,  downward  beams  are  often  seen  at  the  edges  of  auroral 
arcs  (Anderson.  1979),  at  energies  below  the  parallel  potential  drop  within 
"inverted  V's"  (Burch  et  al.  1979;  Lin  and  Hoffman.  1979)  and  in  the  cusp 
( Zanetti  et  al,  1981).  Upward-flowing  and  counter-streaming  electron  beams 
have  recently  been  observed  by  Lin  et  al  (1982)  and  by  Sharp  et  al  (1980). 
Sharp  et  al  attributed  their  S3-3  results  to  fluctuating  double  layers,  and 
reported  that  the  width  in  pitch  angle  of  the  upward  electron  beams  increases 
with  increasing  electron  energy.  We  have  found  that  this  property  is  a 
consequence  of  a  parallel  potential  drop  below  the  spacecraft  and  that  the 
amount  of  spread  is  related  to  the  altitude  of  the  top  of  the  potential  drop. 

Similar  upward  electron  beams  were  reported  at  lower  (1400  km) 
altitudes  by  Johnstone  and  Winningham  (1982)  and  by  Klumpar  and  Heixkila 
(1982).  The  latter  study  showed  that  the  upward  beams  were  imbedded  in  a  much 
larger  region  of  downward  Birkeland  current.  Using  S3-3  data,  Collin  et  al 
(1982)  showed  that  the  occurrence  frequency  of  upward  beams  is  roughly  eight 
times  higher  above  6000  km  than  below  3000  km.  Thus,  one  can  infer  a  sporadic 
low-  altitude  mechanism  and  a  more  stable  mechanism  near  1  R  . 


The  DE-1  spacecraft  is  an  ideal  platform  from  which  to  study  these 
phenomena  for  three  reasons;  (1)  The  spacecraft  moves  slowly  enough  through 
the  acceleration  regions  that  spin  aliasing  is  not  a  problem,  as  it  was  in  the 
Maier  et  al  and  Klumpar  and  Heikkila  studies;  (2)  The  apogee  (~5  R  geocen¬ 
tric)  is  high  enough  that  effects  from  both  the  low-altitude  and  -f  R  accel¬ 
eration  regions  can  be  observed,  as  well  as  any  regions  in  the  1-4  R  altitude 
range;  and  (3)  the  presence  of  upward  electron  acceleration  mechanisms  below 
the  DE-1  altitude  brings  the  cold  ionospheric  electrons  that  carry  the  down¬ 
ward  current  at  low  altitudes  up  into  the  suprathermal  energy  range  where  they 
can  readily  be  detected. 
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Upward  electron  beams  have  been  observed  at  altitudes  above  15,000  km 
on  DE-1 .  These  observations  were  made  primarily  near  the  eguatorward  edge  of 
the  morningside  polar  cusp  (statistically  a  locus  of  downward  current)  and 
sporadically  throughout  the  cusp  region.  Upward  electron  beams  are  seen  all 
through  the  region  of  downward  current,  and  the  magnitude  of  the  currents  car¬ 
ried  by  the  15  eV  to  17  kev  electrons  are  more  than  adequate  to  reproduce  the 
east-west  magnetic  perturbations  measured  concurrently  in  both  the  upward  and 
downward  current  regions.  The  beams,  as  in  other  studies  ( Sharp  et  al,  1980; 
Klumpar  and  Heikkila,  1980;  Collin  et  al,  1982)  ,  are  associated  with  ion 
conics.  The  association  of  these  beams  with  electrostatic  hiss  is  established 
in  another  related  paper  ( Lin  et  al,  1983). 

The  High  Altitude  Plasma  Instrument  (HAPI)  on  Dynamics  Explorer  1 
(DE-1)  made  differential  flux  measurements  of  electrons  and  positive  ions  over 
an  energy/charge  range  of  5  eV/e  to  31  keV/e.  Details  of  the  instrument  and 
its  various  operating  modes  can  be  found  in  Burch  et  al  (1981).  For  the  first 
few  months  after  launch  the  4.7  R£  apogee  of  DE-1  was  just  poleward  of  the 
pre-noon  polar  cusp.  A  common  feature  of  the  HAPI  plasma  data  acquired  during 
this  time  period  is  the  appearance  of  upward- travel ling  electron  beams  with 
energies  from  a  few  tens  of  eV  to  a  few  hundred  eV.  Although  both  upward  and 
downward  electron  beams  are  observed  in  and  near  the  cusp  region,  we  have  con¬ 
sistently  observed  a  region  of  predominantly  upward  beams  that  overlaps  the 
equatorward  boundary  of  the  cusp  electron  population. 

The  systematic  appearance  of  upward  electron  beams  near  the 
equatorward  cusp  boundary,  their  absence  in  the  heart  of  the  cusp,  and 
subsequent  appearance  near  the  cusp  poleward  boundary,  raise  the  question  of 
to  what  extent  the  electrons  measured  on  DE-1  are  charge  carriers  for  dayside 
region  1  and  cusp-region  Birkeland  currents.  This  question  is  answerable  to 
some  degree  by  suitable  correlations  between  the  HAPI  data  and  data  from  the 
DE-1  magnetometer  (Farthing  et  al^,  1981).  The  results  of  the  correlations  for 
a  typical  pass  are  presented  in  the  next  section. 

B.  Electron  Currents  and  Magnetic  Perturbations 

The  parallel  (field-aligned)  current  density  carried  by  a  measured 
electron  distribution  f(v,  a,  d)  is  given  by 

v 
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where  a  is  pitch  angle,  e  is  electron  charge,  v^  and  v^  are  the  velocities 

parallel  and  perpendicular  to  B,  v  .  and  v  are  the  velocity  limits  on  the 

min  max 

measurement,  and  $  is  the  azimuthal  angle  about  B.  since  the  electron  fluxes 

at  energies  below  about  15  eV  are  dominated  by  spacecraft  photoelectrons,  we 

set  vmin  equal  to  the  velocity  of  an  18  eV  electron  (the  next  energy  step 

above  15  eV)  in  the  current  calculations.  In  each  case,  v  was  set  to 

max 

correspond  to  the  maximum  energy  sampled  by  the  HAPI  instrument. 


If  the  spacecraft  potential  is  nonzero,  the  above  calculated 
current  is  an  overestimate,  since  the  measured  electron  velocity  is  larger 
than  the  actual  velocity.  For  the  analyzer  viewing  perpendicular  to  the  spin 
axis  (therefore  also  normal  to  the  spacecraft  surface)  the  velocity  v  in  Eq. 
(4)  should  be  replaced  by  (v2  +  2q<$>1  2,  where  q  =  -e(e)  for  electrons  (ions). 
For  a  potential  4>0  of  roughly  15V,  and  a  total  potential  drop  of  3u  V 
(magnetospheric  and  spacecraft),  we  can  overestimate  the  current  by  no  more 
than  a  factor  of  2.8  and  most  likely  less. 

Calculation  of  the  electron  current  was  performed  as  follows.  A 
differential  current,  dl/dv,  was  first  obtained  for  each  energy  sampled  by 
assuming  a  gyrotropic  distribution  and  applying  a  spline  integration  over  a. 
Next,  the  total  current  was  obtained  by  in  egrating  over  v  with  the  same 
technique.  Ion  current  was  also  computed.  However,  in  all  cases  the  ion 
current  amounted  to  less  than  1%  of  the  electron  current. 

The  results  of  the  current  calculations  for  a  cusp  pass  on  October  b, 
1981,  are  plotted  in  Figure  20,  in  which  upward  currents  are  plotted  as 
positive  values.  Also  plotted  in  Figure  20  arc  the  currents  determined  from 
the  magnetic  perturbations  (dB  )  measured  by  the  DE-1  magnetometer,  assuming 
current  sheets  that  are  greatly  extended  in  longitude  and  oriented 
perpendicular  to  the  local  dipole  magnetic  meridian.  The  agreement  between 
the  computed  and  measured  values  of  the  current  argues  strongly  that 
electrons  with  energies  above  15  eV  are  the  primary  charge  carriers  of  the 
downward  region-1  currents,  the  upward  cusp-region  currents,  and  another 
region  of  downward  current  located  near  the  poleward  boundary  of  the  cusp  in 
the  pre-noon  hours  at  altitudes  near  20,000  km.  These  same  three  current 
systems,  and  corresponding  electron  populations  have  been  observed  in  all  four 
cusp  passes  that  we  have  examined  in  detail  to  date. 

C.  Acceleration  Mechanism 

Detailed  study  of  the  velocity-space  distributions  of  the  observed 
upward  electron  beams  points  to  an  acceleration  region  consisting  of  a 
downward  field-aligned  potential  drop  below  the  spacecraft.  An  example  of  the 
f(v  ,  v  )  distribution  observed  by  the  HAPI  0°  analyzer  on  October  b,  1981,  is 
plotted^in  Figure  21.  To  obtain  optimum  sampling  in  velocity  space.  Figure  21 
contains  log10  f  contours  for  two  successive  spacecrjft  rotations.  The  dots 
in  Figure  8  identify  the  (v  ,  v  )  values  that  were  sampled  by  the  0°  analyzer 
during  the  two  rotations.  To  aicentu.te  the  features  of  the  low  energy  beam. 
Figure  8  contains  data  only  for  energies  up  to  500  eV. 

As  shown  by  Knight  (1973),  conservation  of  energy  and  of  the  first 
adiabatic  invariant  require  that  electrons  originating  in  the  ionosphere 
satisfy  the  inequality  where  V0  is  the  electrostatic  potential  difference 
between  the  ionosphere  aud  the  point  of  observation,  and  B0  and  B  represent 
the  magnetic  induction  at  the  top  of  the  potential-drop  region  and  at  the 
point  of  observation,  respectively.  This  inequality  defines  a  region  within  a 
hyperbola  with  origin  at  (0,0)  and  "semi-axes"  /(2e/m)V0  (intercept  along  the 
Vj|  axis)  and  /(2e/Am)V0,  where  A  =  (B0/B)-l  (asymptotes  v  =  >v  /A1  2).  Shown 
in  Figure  21  is  such  a  hyperbola  for  a  potential  drop  {V 0\  of  31  V  extending 
to  a  maximum  altitude  of  5870  km.  The  interior  of  this  hyperbola,  or  one 
nearly  like  it,  does  in  fact  define  approximately  the  region  occupied  by  the 
upward  electron  beam.  Examination  of  a  number  of  upward  beams  shows  that 
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Figure  21.  Electron  distribution-function  contours  in  the  v]|  ,  vj^  plane 
averaged  over  two  rotations  of  the  DE-1  spacecraft.  Also 
plotted  is  a  hyperbola  defined  by  eq.  (2)  for  a  potential 
drop,  VQ  =  31  V,  with  maximum  altitude  at  5870  km. 
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these  values  for  ve  and  the  altitude  of  the  top  of  the  potential  drop  region 
are  typical.  A  nonzero  spacecraft  potential  $0  will  subtract  from  the 
inferred  accelerating  potential  V0. 

D.  Conclusions 

Upward  electron  beams  are  observed  to  be  a  common  phenomenon  in 
regions  of  downward  Birkeland  current  at  altitudes  above  ~1  R^.  The  electron 
velocity  space  distributions  are  consistent  with  parallel  potential  drops  of  a 
few  tens  of  volts  located  near  1  R£  altitude,  or  if  there  is  a  distributed 
potential  drop  instead  of  a  double  layer,  the  upper  boundary  is  at  about  1  R£ 
and  the  lower  boundary  is  unknown. 

In  the  cusp  region,  electric  current  densities  calculated  from  the 
measured  electron  distribution  functions  over  the  energy  range  18  eV  -  17.5 
keV  are  in  very  good  agreement  with  the  current  inferred  from  the  onboard 
magnetometer  for  both  downward  and  upward  currents.  Uncertainties  in  the 
orientation  of  tne  current  sheet  and  in  the  spacecraft  potential  can  readily 
explain  the  small  differences.  Thus,  while  at  lower  altitudes  downward 
currents  are  typically  carried  by  near-thermal  upgoing  electrons  (with  only 
occasional  upward  beams),  at  altitudes  above  ~1  R£  the  current  is  typically 
carried  by  upward  beams  of  electrons  that  have  been  accelerated  to 
suprathermal  energies  by  potential  drops  of  a  few  tens  of  eV.  Computations  of 
field-aligned  currents  from  hot-plasma  distributions  observed  in  numerous 
orbits  of  DE-1  in  both  the  dayside  and  nightside  auroral  zones  and  in  the 
polar  cap  indicate  that  at  altitudes  near  and  above  10,000  km  the  majority  of 
the  current  carriers  are  electrons  within  the  energy  range  of  the  HAPI 
instrument. 


OBSERVATIONS  OF  A  HEATED  ELECTRON  POPULATION  ASSOCIATED  WITH  THE 
6300  A  SAR  ARC  EMISSION 

Introduction  and  Observations 


It  was  not  until  the  IGY  period  that  the  existence  of  Stable  Auroral 
Red  (SAR)  arc  features  was  first  reported  (Barbier,  1958).  As  implied, 
identification  of  these  subauroral  arcs  using  ground-based  instrumentation 
relies  primarily  on  the  stability  and  monochromatic  nature  of  the  emissions  at 
6300A,  a  result  of  the  1D+3P  transition  of  atomic  oxygen. 

Attention  was  directed  toward  the  ring  current  as  a  potential  source 
of  energy  for  maintenance  of  the  emission  when  early  findings  of  Rees  and 
Akasofu  (1963)  indicated  a  strong  correlation  between  the  presence  of  SAR  arcs 
and  stormtime  enhanced  Dst  indices.  Thermal  conduction  models,  first  proposed 
by  Cole  (1965),  rely  on  the  preferential  excitation  of  atomic  oxygen  to  the  ID 
state  by  electrons  in  the  high  energy  tail  of  a  Maxwellian  distribution.  High 
temperatures  within  the  ionospheric  electron  gas  are  then  maintained  via  con¬ 
duction  of  heat  along  magnetic  field  lines  threading  the  plasmapause-ring  cur¬ 
rent  region.  Several  mechanisms  have  been  suggested  for  the  actual  energy 
transfer.  These  include  the  model  of  Cole  (1965)  which  uses  coulomb  colli¬ 
sions  between  ring  current  protons  and  plasmaspheric  electrons;  that  of 
Cornwall  et  al.  (1971)  which  predicts  the  production  of  ion  cyclotron  wave 
turbulence  within  the  ring  current  and  its  transfer  to  plasmaspheric  elec¬ 
trons;  and  the  more  recent  concept  of  Hasegawa  and  Mima  (1978)  that  introduces 
kinetic  Alfven  waves  to  precipitate  directly  low  energy  electrons  into  the  SAR 
arc  region. 

Recent  observations  by  the  Dynamics  Explorer-2  plasma  experiment  have 
for  the  first  time  provided  measurements  of  the  heated  electrons  responsible 
for  the  production  of  the  6300  SAR  arc  emissions.  In  this  paper  we  will 
present  data  for  a  single  SAR  arc  event  centered  at  04:05  UT  on  day  296  (23 
October)  of  1981.  The  event  was  selected  for  early  analysis  for  two  reasons: 
1)  DE-1  and  DE-2  were  magnetically  conjugate  at  this  time;  2)  the  satellite 
track  was  within  the  window  of  the  Pacific  Northwest  Scanning  Photometer 
(MASP)  network  allowing  detailed  comparisons  between  the  6300A  SAR  emission 
and  the  particle  fluxes.  It  is  in  no  way  a  unique  event. 

On  October  23,  1981,  following  four  days  of  variable  SAR  arc  activity 
as  monitored  by  the  MASP  network,  the  path  of  DE-2  took  it  over  a  distinct  SAR 
arc  during  a  period  and  at  a  location  being  viewed  by  a  MASP  unit  located  at 
Richland,  WA  (46.4°N,  240. 4°E).  Presence  of  the  SAR  arc  was  confirmed  by 
additional  sightings  by  MASP  units  located  near  Eureka,  MT  and  Boulder,  CO. 

Due  to  weather  conditions  and  the  orbital  position  of  DE-2,  we  will  utilize 
data  primarily  from  the  Richland,  WA  site  in  this  paper.  A  detailed  descrip¬ 
tion  of  the  MASP  operation  is  found  in  Kleckner  et  al  (1981). 


Fig.  22  presents  contours  of  6300A  emission  intensity  as  measured 
from  Richland  at  04:00  UT,  near  the  time  of  DE-2  passage.  Standard 
triangulation  procedures  with  the  MASP  unit  located  near  Eureka  indicate  a 
region  of  maximum  emission  at  an  altitude  of  425-450  km.  Compensation  for 
slant  path  viewing  of  the  SAR  arc  feature  and  atmospheric  extinction  have  been 
applied  accordingly.  Also  included  in  Figure  22  is  the  DE-2  track  indicated 
by  the  dashed  line.  The  small  set  of  circular  contours  near  the  center  of  the 


plot  are  artificial,  serving  only  to  locate  the  Richland  site.  The  SAR  arc  is 
clearly  distinguished  by  the  band  of  6300  emission  lying  just  below  L=3. 

Aurora  poleward  of  the  SAR  arc  can  be  seen  beginning  near  L=4. 

DE-1  and  DE-2  were  magnetically  conjugate  during  the  traversal  of  the 
SAR  arc,  giving  a  two  point  simultaneous  plasma  measurement  at  850  and  6000 
km.  The  plasma  experiments  on  board  both  satellites  are  identical  in 
instrumentation  and  response — differing  only  in  configuration.  The  DE-2 
plasma  experiment  is  a  family  of  15  electron  and  15  ion  sensors  mounted  on  a 
single  axis,  magnetically  oriented  scan  platform  allowing  each  detector  to  be 
held  approximately  fixed  in  pitch  angle.  The  DE-1  plasma  instrument  consists 
of  5  electron  and  5  ion  detectors  which  scan  in  pitch  angle  with  satellite 
rotation.  A  complete  description  of  both  instruments  is  found  in  Winningham 
et  al.  (1981)  and  Burch  et  al.  (1981). 

Many  of  the  more  important  plasma  features  of  the  low  energy 
electrons  during  the  SAR  arc  can  be  seen  in  Fig.  23.  The  upper  panel  shows  a 
typical  field-aligned  electron  distribution  function  plotted  in  the  plasma 
rest  frame.  The  electron  population  responsible  for  the  SAR  arc  is  located 
below  8  eV.  The  fact  that  the  distribution  can  be  characterized  by  a  linear 
relationship  in  this  region  suggests  that  a  Maxwellian  fit  to  give  both 
temperature  and  density  may  be  applied.  The  lower  panel  shows  the  complete 
electron  distribution  function  (out  to  50  eV).  The  presence  of  a  field- 
aligned  flow  is  indicated  by  the  offset  of  the  inner  four  contours.  Use  of 
the  offset  (about  275  km/sec  in  this  case)  allow  for  individual  spectra  to  be 
plotted  in  the  plasma  rest  frame.  Calculations  of  both  the  temperature  and 
heat  flux  necessitate  a  transformation  into  the  plasma  rest  frame.  Neglecting 
this  transformation  could  result  in  significant  errors. 

B.  Discussion 

The  identification  of  the  low  energy  electron  population  in  Fig.  23 
as  the  particles  responsible  for  the  6300A  SAR  arc  emission  is  based  primarily 
on  the  correlation  between  the  MASP  photometer  data  and  the  DE-2  particle 
data.  Fig.  24  shows  the  profile  of  the  6300A  emission  along  the  satellite 
track  (solid  line  from  the  Richland,  WA  MASP  unit  and  dashed  line  from  Eureka, 
MT  unit)  in  Fig.  22  along  with  coincident  particle  data.  The  MASP  data 
clearly  show  the  SAR  arc  centered  at  52°A.  The  more  poleward  aurora,  beyond 
56°A,  is  shown  assuming  a  300  km  emission  height,  typical  of  0(1D)  emissions 
in  the  auroral  zone  (dashed  line)  and  at  a  SAR  emission  height  of  420  km 
(solid  line).  The  particle  data  are  shown  in  Fig.  24  as  the  log  of  the 
electron  distribution  function  (7.5°  pitch  angle)  at  6.7  eV  (open  circle), 

8.8  eV  (solid  circle)  and  20.4  eV  (solid  triangle).  In  the  rest  frame  of  the 
plasma  the  lower  two  energies  are  approximately  4.0  and  5.7  eV  respectively. 
Changes  in  the  6.7  eV  values  can  be  interpreted  as  changes  in  density  while 
changes  in  the  6. 7-8. 8  eV  separation  as  changes  in  the  plasma  temperature. 

The  latter,  of  course,  assumes  a  Maxwellian  low  energy  electron  population  and 
as  such  may  not  be  valid  outside  the  general  SAR  arc  region. 

The  correlation  between  the  photometer  and  low  energy  electron  data 
is  unmistakable.  The  distribution  function  at  6.7  eV  tracks  the  6300A 
emission  almost  one  to  one.  The  large  6. 7-8.8  eV  separation  within  the  SAR 
arc  itself  suggests  a  very  cool  electron  population.  Computation  of  the 
electron  temperature  near  52°A  latitude  (taking  into  account  the  plasma  flow 
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Typical  field  aligned  electron  spectrum  taken  near  04:05:20  UT  along 
with  the  corresponding  2-D  phase  space  contour  plot.  The  latter 
distinctly  shows  the  field  aligned  flow  velocity  of  the  low  energy 
electron  population. 


Figure  23. 
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velocity)  gives  values  of  9900  >110*  K.  Corresponding  computations  of  the 
density  using  a  satellite  potential  of  -1.7  ±  0.1  volts  (R.  Heelis,  private 
communication)  gives  a  density  of  250  ±50  cm-*  for  this  particle  population. 

It  should  be  noted  that  there  is  nearly  a  500  km  difference  between 
the  measurement  of  the  particles  and  the  measurement  of  the  6300/  emission 
altitude.  Some  degradation  of  the  particles  in  moving  to  the  emission 
altitude  is  expected,  with  a  single  measurement,  the  amount  of  degradation 
cannot  be  determined  directly  without  a  model.  Without  resorting  to  a 
specific  model,  however,  we  can  make  some  positive  statements  concerning  the 
electron  transport.  Between  days  288  and  324  of  1981  more  than  75  distinct 
particle  populations  with  characteristics  similar  to  those  presented  in  this 
paper  have  been  identified.  When  ground  observations  exist  to  support  the 
satellite  observations  these  electron  populations  have  been  found  to  be 
associated  with  SAR  arcs.  The  satellite  observations  extend  from  above  900  km 
to  400  km  which  indicates  t.*e  ability  of  this  electron  distribution  to  exist 
over  a  considerable  length  of  the  flux  tube. 

The  day  296  SAR  arc  contains  several  features  which  may  yield 
important  information  as  to  the  production  mechanisms  responsible  for  the  low 
energy  electrons.  In  Fig.  24  one  finds  a  decrease  in  the  conjugate 
photoelectron  intensity  within  the  SAR  arc.  Photoelectrons  give  a  measure  of 
the  coupling  between  the  northern  and  soutnern  hemispheres.  Lee  et  al  (1980) 
have  postulated  that  such  decreases  might  be  due  to  enhancements  in  the 
thermal  electrons  within  the  local  flux  tube.  This  shielding  phenomena  has 
been  noted  in  many  of  the  SAR  arc  cases  (both  visual  and  subvisual)  we  have 
investigated. 

There  appears  to  be  no  corresponding  ion  precipitation  accompanying 
the  low  energy  electron  enhancement.  This  is  contrary  to  the  Cornwall  et  al 
(1971)  prediction  of  precipitating  ions  on  the  order  of  10*/cm2-sec.  This 
fact  in  conjunction  with  the  lack  of  a  definite  low  energy  electron  signature 
or  enhancement  during  this  time  period,  might  suggest  a  lower  altitude 
production  for  the  thermal  electrons.  Both  of  these  facts  would  seem  to  point 
to  a  production  mechanism  in  line  with  the  Hasagawa  and  Mima  (1978)  work. 

C.  Summary 

The  low  energy  electrons  observed  by  the  DE-2  plasma  experiment  are 
associated  with  the  6300A  SAR  arc  emission.  The  following  features  have  been 
noted: 

1.  The  heated  electrons  have  a  significant  field  aligned  flow 
velocity. 

2.  There  is  no  corresponding  electron  population  observed  by  the  DE- 
1  plasma  instrument  located  at  6000  km  above  the  SAR  arc. 

3.  The  SAR  arc  occurs  at  the  plasmasheet-plasmasphere  overlap 
boundary,  identified  by  the  drop  out  of  the  plasmasheet  particles 
observed  at  DE-1. 

4.  The  heated  electron  component  appears  fairly  Maxwellian  with 
temperatures  of  about  0.9  eV,  which  is  not  significantly  different 
than  predicted  by  several  theories. 
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5.  The  absence  of  observable  precipitating  ring  current  ions  within 
the  SAR  arc,  which  may  point  to  mechanisms  not  involving  ring  current 
ions. 

6.  The  decrease  in  the  conjugate  photoelectron  intensities, 
indicating  a  partial  shielding  of  one  hemisphere  from  another. 


7.  Modeling  of  the  energy  deposition  and  degradation  of  the  observed 
electrons  should  include  the  possibility  of  a  continual  driving 
source  along  the  total  flux  tube. 


VI.  POLAR  RAIN  OBSERVATIONS 


Winningham  and  Heikkila  (1974)  first  described  polar  rain,  using  data 
from  the  ISIS  1  satellite.  Polar  rain  is  the  most  common  type  of  particle 
precipitation  over  the  polar  caps.  This  broad  and  relatively  structureless 
electron  precipitation  can  often  fill  the  entire  poLar  cap.  The  precipitating 
electrons  typically  have  thermal  energies  on  the  order  of  100  eV  and  are 
isotropic  over  the  downcoming  hemisphere.  The  energy  flux  carried  by  the 
particles  is  of  the  order  of  10-3  to  10-2  erg  cm-2sec-1  (Burke,  1982).  The 
spectral  distribution  of  the  polar  rain  electrons  has  the  same  shape  as  that 
of  cusp  electrons,  but  is  lower  in  intensity,  suggesting  that  polar  rain 
originates  in  the  magnetosheath  and  travels  to  the  polar  ionosphere  via  the 
lobes  of  the  magnetotail. 

Figures  25a  and  25b  show  differential  energy  flux  from  the  LAPI 
instrument  aboard  DE-2  for  a  portion  of  a  polar  pass  on  day  295  (22  October) 
of  19bl.  The  top  panel  shows  data  for  8°  pitch  angle  electrons,  and  the 
bottom  panel  shows  45°  pitch  angle  ion  data.  The  upper  center  panel  shows  90° 
and  0°  pitch  angle  electron  data  from  the  GM  tube.  The  data  are  coded 
according  to  the  color  bars  at  the  right.  Satellite  ephemeris  is  shown  at  the 
bottom.  In  this  particular  pass,  DE-2  passed  through  the  northern  polar  cusp 
at  approximately  1306-1307  UT  and  into  the  polar  cap.  The  polar  rain  can  be 
seen  clearly  as  a  bright  band  in  the  electron  panel  over  the  energy  range  ~60- 
600  eV.  The  electron  flux  intensity  remained  fairly  constant  in  intensity 
until  just  after  1309  UT  when  it  abruptly  decreased.  Notice  the  lack  of  ion 
fluxes  over  the  polar  cap. 

Figure  26  shows  line  plots  of  the  average  energy,  energy  flux,  and 
density  of  the  downcoming  electrons  for  1306-1311  UT.  These  parameters  are 
obtained  by  integrating  distribution  functions  over  energy  from  5  eV  to  20 
keV.  During  the  entire  interval,  the  average  electron  energy  and  density 
remained  roughly  constant.  The  average  energy  was  about  100  eV  and  the 
density  of  downcoming  electrons  was  in  the  range  of  1-5  cm-3.  These  values 
of  energy  and  density  are  representative  of  the  polar  rain.  After  exit  from 
the  cusp  at  about  1307  UT,  the  energy  flux  increased  very  slightly  toward  the 
pole  until  gust  after  1309  UT  when  it  decreased  by  approximately  50%. 

Figure  27  shows  a  DE-1  high-altitude  polar  pass  approximately  three 
hours  after  the  DE-2  pass  shown  in  Figure  25.  During  the  interval  lbl7-1619 
UT,  DE-1  traversed  the  magnetic  local  time  and  invariant  latitude 
corresponding  to  those  of  DE-2  just  poleward  of  the  cusp  in  Figure  25a.  The 
format  is  the  same  as  for  Figure  25,  except  for  the  absence  of  GM  tube  data. 
The  polar  rain  may  be  seen  as  a  band  in  the  spectrogram  over  the  energy  range 
~ 100-b0o  eV.  Bracketing  the  polar  rain  are  counterstreaming  energetic  (E>600 
eV )  electron  beams  at  0°  and  180°  pitch  angle  and  upward-moving  low-energy 
( E< 100  eV)  electron  beams.  The  source  of  the  energetic  counter streaming 
electrons  is  open  to  question  and  presently  under  study  by  Gurgiolo  and  Burch 
(1982).  Two  possibilities  for  the  source  would  be  reflection  by  a  potential 
barrier  at  the  magnetopause  (Foster  and  Burrows,  1976)  or  direct  access  of 
solar  electrons.  Upstreaming  ion  conics  may  be  seen  in  the  lower  panel. 

Comparing  Figure  25  and  Figure  27,  it  is  interesting  to  note  that  the 
polar  rain  is  of  the  same  spectral  nature  at  both  altitudes,  suggesting  the 


absence  of  any  significant  acceleration  between  the  two  satellite  locations. 

To  demonstrate  farther  the  similarity  between  the  electrons  at  the  two 
satellites,  Figure  28  shows  number  density,  energy  flux,  and  average  energy  at 
DE-1  for  the  period  1614-1619  UT.  The  average  energy  was  ~10Q  eV,  density  was 
-1  cm-*  and  the  energy  flux  was  0.8  ergs  cm-2s-1.  The  corresponding  values  at 
lower  altitudes  were  -100  eV,  -2  cm-*,  and  -0.8  ergs  cm-2s-1. 

A  detailed  look  at  the  near-field-aligned  electron  populations  at  the 
two  satellite  locations  is  provided  in  Figure  29.  Shown  are  three  data 
sets  -  two  single  energy  sweeps  from  LAPI,  one  just  at  the  poleward  edge  of 
the  cusp  and  one  approximately  10°  poleward  of  the  cusp,  and  a  single  energy 
sweep  from  HAPI  at  a  location  which  maps  to  near  the  location  of  the  first 
LAPI  spectrum.  From  Figure  27,  we  first  can  see  that  above  photoelectron 
energies  the  LAPI  population  can  be  represented  by  a  Maxwellian,  allowing  the 
temperatures  to  be  estimated.  Near  the  cusp  a  least  squares  fit  gives  a 
temperature  of  92  eV,  while  farther  poleward  the  temperature  is  99  eV.  The 
HAPI  spectrum  appears  well-fitted  by  two  Maxwellians  above  photoelectron 
energies.  The  low-energy  population  has  a  least-squares  fit  temperature  of 
65  eV,  while  the  higher  energy  population  has  a  temperature  of  141  eV. 
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FIGURE  29.  Electron  energy  spectra  from  DE-1  and  DE-2  in  the  cusp  region 
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VII.  ACCELERATED  POLAR  WIND 


The  polar  wind  refers  to  the  continual  escape  of  ionospheric  plasma  from 
the  polar  ionosphere  along  open  magnetic  field  lines.  This  escape  of 
particles  has  consistently  been  observed  to  lead  to  a  depletion  of  light 
ionospheric  ions  (ref.  7)  and  electrons  in  the  upper  polar  ionosphere.  Until 
recently,  observations  of  the  polar  wind  in  the  literature  were  limited  to 
the  low-altitude  observations  of  Huffman  an^  Dodson  (1980).  They  reported  a 
continual  upward  flow  of  H  and  He  over  the  polar  caps  in  the  range  1-5 
km/sec  and  1-3  km/sec,  respectively.  The  low  energy  of  these  particles 
precluded  most  instruments  from  directly  obtaining  a  distribution  function. 
However,  as  pointed  out  by  Gurgiolo  and  Burch  (1982),  polar-wind  models 
predict  that  the  particle  velocity  snould  increase  with  altitude  (Schunk  and 
Watkins,  1982).  Thus,  the  high  altitude  satellite  DE-1  is  an  ideal  platform 
from  which  to  study  the  polar  wind. 

Plasma  data  from  DE-1  polar  passes  indicate  that  ions  with  peak 
differential  energy  fluxes  in  the  5  to  100  eV  range  are  continually  flowing 
out  of  the  dayside  cusp  and  polar  cap.  The  flows  have  both  a  field-aligned 
and  a  conic  component.  The  field-aligned  component  is  unmistakably  the  polar 
wind.  Gurgiolo  and  Burch  (ref.  5)  concluded  that  the  conics  observed  in 
conjunction  with  the  polar  wind  are  polar-wind  ions  that  have  been 
perpendicularly  heated. 

The  polar-wind  observations  to  be  discussed  were  made  during  a  polar 
pass  on  day  272  (September  29)  of  1981.  Figure  30  consists  of  3  sets  of 
particle  spectrograms  on  the  same  format  as  Figure  3.  Each  spectrogram 
displays  the  differential  energy  flux  for  the  particle  detectors  which  lie 
closest  to  the  plane  containing  the  local  magnetic  field  vector. 

Figure  30  shows  a  continuous  band  of  upward  flowing  ions  (near  180° 
pitch  angle).  The  satellite  was  inside  the  cusp  during  the  interval  covered  in 
Figure  30a.  At  1411  UT,  the  satellite  passed  through  the  poleward  cusp  bound¬ 
ary.  The  upward  ions  showed  a  gradual  increase  in  energy  throughout  and 
slightly  poleward  of  the  cusp,  then  they  begin  to  decrease  steadily,  the  peak 
in  the  differential  energy  flux  eventually  dropping  below  the  lowest  energy 
channel  of  HAPI.  Within  the  cusp,  the  conic  and  field-aligned  ions  appear  as 
"tuning  forks"  in  the  spectrograms.  In  the  polar  cap  the  conics  become  less 
and  less  apparent,  as  the  energy  and  angular  separation  of  the  conic  and 
field-  aligned  ions  decrease.  In  the  polar  cap,  there  also  is  a  shift  in  the 
field-  aligned  ions  from  180°  pitch  angle  to  lb5°,  indicating  an  antisunward 
flow. 


Figure  31  shows  a  detailed  view  of  the  field  aligned  ion  population.  Two 
data  sets  are  shown,  each  of  which  was  averaged  over  four  satellite 
revolutions  to  improve  counting  statistics.  Each  of  the  distributions  was 
transformed  to  the  rest  frame  of  the  plasma  prior  to  averaging  (except  for  S/C 
charging  which  was  about  20V).  The  two  data  sets  were  taken  beginning  1434:05 
UT  (solid  circles)  and  1434:54  UT  (open  circles).  In  Figure  31  Maxwellian 
distributions  are  straight  lines,  with  the  slope  being  proportional  to  the 
plasma  temperature  and  the  intercept  with  the  distribution  axis  being 
proportional  to  the  plasma  density. 


From  Figure  31  we  can  see  that  the  field-aligned  ion  component  observed 
in  the  polar  regions  is  comprised  of  two  ion  populations--a  low  energy 
component  with  temperature  below  0.5  eV  and  a  high  energy  tail  with 
temperature  above  1.5  eV.  Both  ion  populations  would  appear  to  be  well 
represented  by  a  Maxwellian  distribution,  allowing  the  estimates  of  the 
temperatures  to  be  made.  A  least  squares  fit  to  the  low-energy  population  of 
Figure  7  gives  a  temperature  of  0.29±0.16  eV  ( 320Q°K±1800K ) .  The  corresponding 
temperature  of  the  high  energy  plasma  population  in  Figure  7  is  found  to  be 
2.7±0.7  eV.  The  peak  at  about  1  eV  may  be  due  to  an  ion  heavier  than  H  . 

There  is  little  doubt  that  the  lower  energy,  field-aligned  ions 
constitute  the  "classical"  polar  wind.  The  high  energy  tail,  however,  is  not 
a  feature  predicted  by  polar-wind  models.  It  is  likely  that  this  hotter 
plasma  is  the  result  of  a  perpendicular  heating  of  polar-wind  ions.  Assuming 
that  all  ions  in  the  polar  wind  are  of  equal  temperature,  Figure  7  gives 
heating  on  the  order  of  4  to  10. 

We  envision  a  scenario  as  shown  in  Figure  32  occurring  during  the  polar 
wind  escape  along  magnetic  field  lines.  At  low  altituaes  the  generation  of 
the  polar  wind  occurs  and  ionospheric  ions  begin  to  escape  along  the  open 
polar  cap  and  cusp  field  lines  (Figure  32a).  The  presence  of  ion 
perpendicular  heating  along  the  open  field  lines  produces  the  characteristic 
conic  signature  (Figure  32b).  The  distribution  above  the  heating  region  can 
be  considered  a  two-component  plasma  consisting  of  unheated  and  heated  polar- 
wind  particles.  The  unheated  polar  wind  is  still  field  aligned,  while  the 
heated  ions  have  large  velocities  in  the  direction  perpendicular  to  B.  As  the 
distribution  travels  upward,  the  perpendicular  arms  of  the  distribution  will 
begin  to  collapse  toward  the  magnetic  field  direction  in  accordance  with  the 
first  adiabatic  invariant,  and  a  field-aligned  high  energy  tail  to  the  polar 
wind  should  develop  (Figure  32c).  Figure  32d  shows  contours  of  a  typical 
conic/polar-wind  population  measured  in  the  cusp  on  day  272. 

Ion  conics  have  been  reported  by  numerous  people  (Ungstrup  et  al.,  1979; 
Gorney  et  al.,  1981).  Theories  as  to  the  production  of  the  conic  ion 
distribution  have  favored  a  perpendicular  neating  by  electrostatic  ion- 
cyclotron  waves  (Dusenberry  and  Lyons,  1981;  Lysak  et  al.,  1980),  although 
recently  (Chang  and  Coppi,  1981)  it  has  been  proposed  that  a  heating  by  lower 
hybrid  waves  may  also  be  a  viable  conic  generation  mechanism.  The  problem  is, 
however,  that  most  of  these  theories  are  applicable  only  at  low  altitudes  in 
the  auroral  zone.  Adaptations  of  the  theories  to  the  environment  of  the 
polar  cap  must  be  made  to  explain  these  observations. 

Using  the  observed  conic  pitch  angle  distribution,  the  altitude  of  the 
observation  and  an  assumed  value  for  the  initial  conic  pitch  angle,  it  is  pos¬ 
sible  through  the  first  adiabatic  invariant  to  estimate  tne  altitude  of  the 
conic  generation.  In  general  the  initial  conic  pitch  angle  for  such  computa¬ 
tions  is  assumed  to  be  90°;  however,  if  the  source  plasma  already  has  a 
significant  Vo,  as  the  polar  wind  does,  then  the  initial  conic  pitch  angle 
is  expected  to  shift  to  lower  values.  Using  a  dipolar  magnetic  field 
approximation  and  estimating  the  values  of  the  initial  conic  pitch  angle 
through  energy  conservation  arguments  (which  assume  that  Vo  does  not  change 
between  the  observation  height  and  the  conic  height)  we  estimate  the  heating 
to  occur  at  an  altitude  of  about  12,000  km--constant  throughout  both  the  cusp 
and  polar  cap.  By  using  an  initial  conic  pitch  angle  of  90°  we  can  place  a 
lower  limit  to  the  heating  altitude  at  8000  km. 
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FIGURE  32.  Evolution  of  polar  wind.  Shown  are  isocontours  of  the  distribution 
function  in  phase  space 
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VIII.  OBSERVATIONS  OF  ION  AND  ELECTRON  ACCELERATION  EVENTS 
PRODUCED  BY  PARALLEL  ELECTRIC  FIELDS 

The  electrical  coupling  between  the  high-altitude  and  low-altitude 
regimes  of  the  polar  cap  is  an  important  aspect  of  magnetosphere-ionosphere 
interactions.  The  unique  placement  of  the  Dynamics  Explorer  orbits  with 
respect  to  each  other  has  afforded  an  excellent  opportunity  to  investigate 
this  coupling  between  the  two  regimes.  Frequently,  DE-1  and  DE-2  observe 
signatures  of  field-aligned  acceleration  of  ions  and  electrons  above  the  polar 
cap.  In  this  section,  we  briefly  present  observations  made  on  October  17, 

1981  (day  81290)  above  a  "theta  aurora"  signature. 

Figure  33  shows  a  spectrogram  of  180°  ions  and  0°  pitch  angle  electrons 
from  MAPI.  The  satellite  crossed  the  polar-cap  arc  field  lines  during  the 
interval  1630-1650  UT.  Two  intervals  of  intense  electron  fluxes  below  1  keV 
were  seen:  -1631:30-1633:15  UT  and  -1642-1649  UT.  Figure  34  shows  contour 
plots  of  the  log  of  the  distribution  functions  of  ions  (Figure  34a)  and  elec¬ 
trons  (Figure  34b)  for  the  interval  1632:45-1632:57  UT.  The  ions  showed  a 
strong  upgoing  beam  at  approximately  luO  eV.  The  electron  fluxes  were  iso¬ 
tropic  except  for  a  loss  cone.  Thus,  the  combined  electron  and  ion  observa¬ 
tions  indicated  that  an  approximately  100  eV  potential  drop  lay  below  the 
satellite  altitude. 

Figure  35  shows  contours  in  the  same  format  as  Figure  34,  for  the 
interval  1646:30  to  1646:36  UT.  By  this  time  the  ions  had  become  isotropic, 
while  the  electrons  showed  a  downcoming  beam  of  100  eV.  Thus,  the 
acceleration  region  appears  now  to  be  above  the  satellite  altitude  (17313  km) 
while  remaining  at  -100  eV. 

The  plasma  characteristics  described  above  are  very  similar  to  those 
observed  by  DE-1  in  the  high  altitude  (>15,000  km)  nightside  auroral  zone, 
except  that  the  energy  of  the  polar  cap  acceleration  events  is  lower. 
Correlation  studies  by  Hardy  et  al.  (1982)  have  already  noted  that  the  plasma 
signatures  of  polar  cap  acceleration  events  at  low  altitudes  are  similar  to 
those  of  evening  discrete  arcs.  The  parallel  electric  field  responsible  for 
the  polar-cap  acceleration  could  be  produced  by  the  same  process  as  the 
evening  auroral  arcs. 

Satellite  observations  of  convection  and  electron  precipitation  at  low 
altitudes  (<1000  km)  have  indicated  that  convective  electric  fields  point 
toward  the  region  of  polar-cap  acceleration  (V'E<0)  (Burch  et  al.,  1979).  The 
polar  cap  acceleration  region  is  therefore  a  region  of  negative  space  charge. 
Burke  et  al.  (1982)  used  the  simultaneous  electric  field,  magnetic  field  and 
electron  flux  measurements  of  the  S3-2  satellite  to  demonstrate  that  upward 
Birkeland  currents  were  embedded  in  regions  of  polar  cap  acceleration.  These 
observations  suggest  that  field-aligned  potentials  develop  as  a  result  of 
imperfect  mapping  of  magnetospheric  convective  electric  fields  to  the  iono¬ 
sphere  (Hardy  et  al.,  1982;  Burke  et  al.,  1982). 
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